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ABSTRACT 

Two-dimensional numerical simulations of the formation of cold condensations 
in a vertical current sheet have been performed using the radiative, resistive 
MHD equations with line-tied boundary conditions at one end of the sheet. 
Prominence -like condensations are observed to appear above and below an X-line 
produced by the onset of the tearing-mode instability. Cooling in the sheet is 
initiated by Ohmic decay, with the densest condensations occurring in the region 
downstream of a fast-mode shock. This shock, which is due to the line-tied 
boundary conditions, terminates one of the two supermagnetosonic reconnection 
jets that develop when the tearing is fully developed. This paper emphasizes the 
condensation properties of shock waves, which may trigger or considerably 
enhance the conditions for thermal condensations. 

I. 2D NUMERICAL SIMULATION 

The initial condition is a vertical current sheet in both mechanical and 
thermal equilibrium, and the half-width and the height of the sheet are 0.15 and 
4.0 in units normalized to the horizontal size of the box. The sheet magnetic 
Reynolds number is 120; and the ratios of radiative, diffusive and tearing time 
scales to the Alfv^nic one are 20, 120, and 11, respectively. A symmetry 
condition is used in the center of the sheet, free-floating boundary conditions 
are fixed at the top and the right edges of the box, and line-tying conditions 
are used at the base. The initial plasma g outside the current sheet is 0.1 
(Forbes and Priest 1983, Malherbe et al. 1984, Forbes and Malherbe 1986a, b) . 

We now present an order of magnitude model to investigate the thermal 
effects of hydrodynamic shock waves and explain their condensing properties in 
the solar corona. 

II. A SHOCK CONDENSATION MECHANISM FOR PROMINENCES 

Consider a hot coronal equilibrium (subscript o quantities) described by: 
°-hp o - P^8(T 0 ), (1) 

where the wave heating term hp^ balances the radiative loss term 
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3u/3t = - u [t* 1 - t" 1 (1-0)] = -U/T, 

C R 

2 5/2 

where x^ = (P q L C ) A Q T ) and x^ = (T q C )/[p Q(T q )] are the conductive and 
radiative time-sca¥es .° 8he hot equilibrium's s?able when 

“1 -1 x -1 

t = x - (1- a) T 
c R 

With „no conduction (t = 00 ) this relation is always satisfied for ot > 1 or T <8 
x 10 H K. So a cold equilibrium is always stable and a hot one always unstab?e. 
With conduction, the equilibrium is stable when x /x > 1-a and this is always 
the case when ot > 1 (cold equilibrium) . It may a?so C be the case when a < 1 (hot 
equilibrium) , if the following inequality is satisfied: 

L < L = [k T 7/2 ] lj [p^Q(T Ml-a)]’ 1 * . 
c o o o o 

If L > L , the hot equilibrium becomes unstable, 
c 

Suppose now that a fast-mode MHD shock occurs in the reconnection jet, and 
let us examine the quantities x , x , x, and L , upstream (supersonic flow 
denoted by subscript u) and downstream (subsonic flow denoted by subscript d) of 
the shock. Because the magnetic field in the jet has almost been completely 
annihilated we can use the hydrodynamic jump relations. For a monatomic gas 
with y = 5/3, these give 

1 ; e a /f u ; 4- 

Defining the compression factor across the shock as x = p^/p^, We ^ ave 

V u /V d = X ' P d /P u = < 4 X-D/(4-X) ' T u /T d = X(4-X)/(4X-1) • 

From (4) , we obtain 

T Ru A rd = x2 " a [(4 "* )/(4 “I)] 1 " 0 

T cu /T =d- X ~ 7/2 (<4X-l)/(4-X>] 5/2 

_ /T ll/4-a/2 ... . ... . . 7/4 -a/2 , , . 

L cu /L cd = X [(4-x)/(4x-l)] , (a < 1) 

5 

At high temperatures (T > 3 x 10 K, a < 0) and in the absence of 
conduction, the shock decrease^ the radiative losses (x r /x^ < 1) • However, at 
lower temperatures (T < 3 x 10 K, a 0) , it increases \:he radiation losses, and 
the cooling time becomes faster. 

When the temperature is high (T > 3 x 10 5 K) , the shock makes the triggering 

of a thermal instability more difficult in the presence of conduction (L /L ) . 

But, when T < 2 x 10 K(a = 0.9) , it is possible to have a thermally sta£¥e C 

upstream (L < L ) and an unstable downstream (L > L ) . This is the case when 

L /L _ > 1 and when L , < L < L . As one can see 5rom Figure 3 this 
cu cd cd cu 

corresponds to a shock strength X smaller than 2 when a = 0.9. Hence, under 
certain conditions (e.g. transition zone like temperatures) the shock may trigger 
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